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The magnetocentrifugal disk wind mechanism is the leading candidate for producing the large-scale, 
bipolar jets commonly seen in protostellar systems. I present a detailed formulation of a global, radially 
self-similar model for a non-ideal disk that launches a magnetocentrifugal wind. This formulation 
generalizes the conductivity tensor formalism previously used in radially localized disk models. The 
model involves matching a solution of the equations of non-ideal MHD describing matter in the disk 
to a solution of the equations of ideal MHD describing a "cold" wind. The disk solution must pass 
smoothly through the sonic point, the wind solution must pass smoothly through the Alfven point, 
and the two solutions must match at the disk/wind interface. This model includes for the first time a 
self-consistent treatment of the evolution of magnetic flux threading the disk, which can change on the 
disk accretion timescale. The formulation presented here also allows a realistic conductivity profile 
for the disk to be used in a global disk/wind model for the first time. The physical constraints on the 
model solutions fix the distribution of the magnetic field threading the disk, the midplane accretion 
speed, and the midplane migration speed of flux surfaces. I present a representative solution that 
corresponds to a disk in the ambipolar conductivity regime with a nominal neutral-matter-magnetic- 
field coupling parameter that is constant along field lines, matched to a wind solution. I conclude 
with a brief discussion of the importance of self-similar disk/ wind models in studying global processes 
such as dust evolution in protostellar systems. 

Subject headings: accretion, accretion disks — ISM: jets and outflows — magnetohydrodynamics — stars: 
formation — stars: protostars — stars: winds, outflows 



1. INTRODUCTION 

Low-mass star formation typically involves the collapse 
of a dense core inside a molecular cloud to form a con- 
densed central object and an extended disk. The pro- 
tostar gains much of its mass by accretion of disk ma- 
terial, which can only proceed if there is some mecha- 
nism for removing angular momentum from the accret- 
ing material. Protostars that are actively accreting disk 
material frequently exhibit associated bipolar outflows 
in the form of coUimated jets oriented perpendicular to 
the disk, stretching out over distances much larger than 
the radial scale of the disk (e.g., Bally et al. 2007). This 
association between accretion disks and large-scale, col- 
limated jets is also seen in other systems, including X- 
ray binaries and active galactic nuclei. The strongest 
candidate mechanism for launching jets in these systems 
is described by the m agnetocentrifugal disk wind model 
(jBlandford fc P avnc 1982. hereafter BP82), which relies 
on the acceleration of a fraction of disk material out along 
a large-scale, ordered magnetic field threading the disk 
at a sufficiently small inclination angle. This magnetic 
field could be the interstellar field in the parent molec- 
ular cloud that is dragged in during core collapse, or it 
could be generated by dynamo activity in the disk or 
the star. A natural consequence of the disk wind model 
is that the wind efficiently removes angular momentum 
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from the disk material, which could explain the strong 
inferred link betw een accretion and outflow phenomena 
in protostars (e.g.. iKonigl &: Salm cron"2011f). 

Several previous semianalytic studies have matched so- 
lutions of the equations of non-ideal MHD that describe 
a weakly ionized, magnetized accretion disk, and solu- 
tions of the equations of cold, ideal MHD for a poly- 
tropic fluid that describe a BP82-type magnetocentrifu- 
gal wind that removes all the disk's excess angular mo- 
mentum. These models f ocused on a radially loc alized 
portion of th e disk fe.g.. iWardle fc Konigll [l993l here- 
after WK93: ISalmeronet all 120071 hereafter SKW07; 
iSalmcron elal| |201,1', hereafter SK Wll) or on a global 
disk mo del Xe.gTlKonigl 1989; 119961 hereafter L96; 
iFerrcira 1997, hereafter F97) based on the same assump- 
tion of self-similarity in spherical radius that underlies 
the BP82 wind model. Another approach is presented in 
iCampbelll (|2003l . I2005D . where the thinness of the disk is 
used to justify solutions for the disk structure that are 
separable in the radial and vertical coordinates; these 
solutions are matched to an isothermal wind. A key pa- 
rameter in determining whether a disk can launch a mag- 
netocentrifugal wind is the field line inclination angle at 
the disk surface _Br,s/5z,sj where _Br,s and Bz^s are the 
radial and vertical field components at the disk surface 
(subscript 's'). This inclination angle is determined by 
the fiux distribution along the disk surface and can only 
be cal culated in a self-consist ent fashion in a global treat- 
ment (jOgilvie fc Livioll200lL hereafter OLOl). A second 
key parameter describing the interaction of the disk and 
wind is -B0,s, the aziniuthal field strength at the disk sur- 
face, which sets the torque exerted by the wind on the 
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disk (BP82). As with the surface field inchnation angle, 
the value of B^^s is determined by processes acting out- 
side the disk; specifically, the value of iJ^^s is determined 
by the constraint that the wind solution pass smoothly 
through the Alfven point. 

Accretion fiows in a weakly ionized disk will drag in 
magnetic flux, opposing the natural tendency of the mag- 
netic flux to diffuse outwards. Studies of accretion disks 
that shed angular momentum due to MHD turbulence 
developed from the magnetorotational instability (MRI, 
iBalbus fc~H awlcv 1991) suggest that the magnetic diffu- 
sivity and the efi'ective viscosity of such disks are compa - 
rable (|Guan fc Gammiell2009l: iLesur fc Longarett|l2009D 
agreeing with previous work bv iLubow et al.l (|1994aD 
which suggested that the magnetic diffusivity of MRI- 
active accretion disks may be too large to allow efficient 
inward dragging of magnetic flux. This problem can be 
overcome in an MRI-active disk threaded by a large-scale 
fleld that removes angular moment um from surface layers 
in which the MRI i s supp ressed (jRothstein fc Lovelacd 
l2008HLovelace et al.|[2009l) . More generally, this problem 
disappears if the angular momentum transport mech- 
anism is decoupled from the mechanism that provides 
the magnetic diffusivity so that the diffusivity is much 
smaller than the effective viscosity of the disk. This oc- 
curs naturally in a disk from which angular momentum 
is removed through a magnetocentrifugal wind, allowing 
magnetic flux to be dragged in and distributed over a 
wind-driving disk. As noted by OLOl, the evolution of 
the magnetic flux through the disk determines whether 
magnetic conflgurations capable of launching outflows 
can be maintained over the accretion timescale. Previous 
studies of wind-driving disks have not treated the effects 
of magnetic flux evolution in a self-consistent manner: 
in a local model, OLOl calculated the rate of flux migra- 
tion given Brs and B^,^, whic h were free p arameters. 
The local mo'dels of iCampbelll (|2003l [20051 ) expressed 
Bcft^s and the surface field inclination angle in terms of 
unconstrained parameters that describe the structure of 
the disk and the processes operating within it. WK93 
also employed a local model in which they determined 
i?0^s self-consistently using the Alfven point constraint, 
but left the surface field inclination angle unconstrained, 
describing the magnetic fiux migration speed using the 
free parameter €b- The global models of L96 and F97 
also solved for B^^g using the Alfven point constraint 
but left the surface field inclination angle and magnetic 
flux migration unconstrained, essentially setting cb = 0. 
For an overview of recent work on semianalytic models 
and complementary n umerical simulations of prot ostellar 
jets, see the article bv lKonigl fc SalmeronI (|2011[ ). 

I present here a new, semianalytic, radially self-similar, 
global model that matches disk solutions with BP82-type 
wind solutions, adopting an approach similar to that of 
L96 but including for the first time a fully self-consistent 
calculation of the effects of magnetic flux migration and 
the magnetic fleld structure above the disk: cb, Br.s, 
B^^s and Br^s/Bz,s are no longer free parameters but are 
determined by physical constraints imposed at the disk 
surface and at critical points of the outflow. In addi- 
tion, the equations describing the disk are formulated in 
terms of the conductivity tensor, allowing for fully gen- 
eral non-ideal disk models. The plan of the paper is as 
follows. In Section [2] I review the equations describing a 



(nearly) steady-state, axisymmetric, weakly ionized disk 
threaded by a magnetic field, and I show how the sys- 
tem reduces to a set of ordinary differential equations 
under the radial self-similarity assumption. In Section [3] 
I discuss constraints from the self-consistent treatment 
of the field-line inclination at the top of the disk and 
from matching disk and wind solutions. In Section 2] I 
present a representative matched disk/wind solution, and 
in Section [5] I summarize the study and outline plans to 
improve on the current model. 

2. FORMULATION 

I consider a weakly ionized, geometrically thin, ver- 
tically isothermal, axisymmetric disk in near-Keplerian 
rotation around a protostar of mass M that is threaded 
by an ordered, large-scale magnetic field. I examine 
the general, non-ideal MHD case fSection l2.1l) . express- 
ing the matter-field interaction using the conductivity 
tensor formalism, and seek solutions that are steady 
on timescales shorter than the accretion time (Sec- 
tion [2TOJ. I follow L96 and F97 in making the sim- 
plifying assumption of radial self-similarity to produce 
a global disk/wind model (Section I2.1.2p . and I derive 
the midplane values of various physical quantities (Sec- 
tion d^^]). 

In Section 12.21 1 specialize to the case of the ambipolar 
conductivity regime. I once again derive the relevant self- 
similar disk equations (Section I2.2.2[) and the midplane 
values for the variables of integration (Section I2.2.3P . 

2.1. General Non-Ideal Case 

The interiors of protostellar disks have a low ionization 
fraction, necessitating the use of non-ideal MHD to de- 
scribe the behavior of the disk material in the presence 
of a large-scale magnetic field that is assumed to thread 
the disk. The non-ideal matter-field interaction can be 
described using a multi-fluid approach or by introducin g 
the conductivity tensor fe.g.. IKonigl fc SalmeronI [20 111) , 
which encapsulates the response of the dominant neutral 
fluid to the presence of the magnetic fleld. The effects 
of the magnetic fleld are communicated to the neutrals 
by the development of charge-neutral drift velocities due 
to the Lorentz forces on charged species, and the conse- 
quent systematic charge-neutral drag forces. Under the 
conductivity tensor formalism, this is expressed in the 
generalized version of Ohm's law, 

J = cr ■ Ec ~ (7o-E'c|| + fH-B ^ Ec± + (TpEc± , (1) 

where Ec is the electric field in a frame comoving with 
the dominant neutral fluid, E^w and Ec± are the compo- 
nents of Ec parallel and perpendicular to the magnetic 
field, respectively, and cris the conductivity tensor. The 
components of crare the Ohm conductivity ctq, the Hall 
conductivity ith , and the Pedersen conductivity tip . 
For a given charged species (denoted by subscript 'j'), 

the product of the gyrofrequency (^^-j^^^^ and the 

neutral-charge momentum exchange time 1/ (7jp) gives 
a measure of the coupling of the charged species to the 
magnetic field in a sea of neutrals. In these expressions, 
Zj, nrij, and 7j are the signed charge, mass, and charge- 
neutral coUisional drag coefficient of species j, and p is 
the neutral mass density. The signed value of the ratio of 
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the gyrofrequency and neutral collision frequency is the 
Hall parameter, 



mjC -ijp 



(2) 



The definition of the Hall parameter uses the signed value 
of the magnetic field strength, B = \B\sgn{Bz}, where 
Bz is the z component of the magnetic field in cylin- 
drical coordinates {r, (j), z}, so that the Hall conductivity 
retains its dependence on the magnetic field polarity (see 
Equation Large values \/3\j ^ 1 imply that charged 
species j is tightly coupled to the magnetic field, whereas 
small values \/3j \ <C 1 imply that species j is weakly cou- 
pled to the field. 

The conductivity tensor components can be expressed 
in terms of the Hall parameters as 



and 



CTH 



X ^ II 



nj\Zj\ 



ec X ^ 



1+/? 



(3) 
(4) 

(5) 



Note that since B is the signed value of the magnetic 
field strength, the Hall conductivity can be positive or 
negative. The factor of B in /3j cancels the pre-factor of 
in the expressions for the Ohm and Pedersen con- 
ductivities, which are positive regardless of the magnetic 
field orientation. 

The coupling of the neutral matter to the magnetic 
field is described by the Elsasser number A, given by 



A 



(6) 



where wa = 



/Anp is the Alfvcn speed, fix is the 



Keplerian angular frequency, 7y_L = c /47r(T_L, and tr^ = 
(0-2 + crlf/^ (e.g., IKonigl fc SalmeronI [ml . As with 
the Hall parameter, large values A ^ 1 imply that the 
neutral matter is strongly coupled to the magnetic field, 
whereas small values A ^ 1 imply that the neutral mat- 
ter is only weakly coupled to the magnetic field. Launch- 
ing a magnetocentrifugal wind typically requires strong 
neutral-field coupling (A > 1). 

2.1.1. Governing Equations 

The governing equations for the disk are the continuity 
of mass, 

V-(pV)=0, (7) 
and the conservation of momentum. 



dV 1 
^ + (F-V)F + -VP 

at p 



J X B 

cp 



0. 



(8) 



In these equations. Vis the neutral gas velocity and P = 
yOCg is its pressure, where Cg is the isothermal speed of 
sound. The gravitational potential of the central object 
is given by $ = -GM/R = -V^{R), where R = {r"^ + 



z^y^^ is the spherical radius and Vk is the Keplerian 
rotation speed. 

Within the disk, the ionization fraction is very low, 
which implies that the equations of momentum conserva- 
tion for the charged species are dominated by the charge- 
neutral drag force and the Lorentz force (cf. Equations 
2.15-2.17 in WK93). We can thus replace the charge- 
neutral drag force with the Lorentz force; this substitu- 
tion has already been made in the conservation of mo- 
mentum for the neutral disk material (Equation (|8])). 

The current density satisfies Ampere's law (neglecting 
the displacement current) 



47r 



V X S, 



(9) 



and Ohm's law (Equation ([T])), with the comoving elec- 
tric field Ec related to the inertial-frame electric field 
Ehy 

Ec = E + ^^. (10) 
c 

The magnetic field satisfies the solenoidal condition 



(11) 



and the electric field £Jsatisfies the induction equation. 



dB 

'dt 



= -cV X E , 



or, in terms of the comoving electric field, 



-^Vx(VxB) 



cV X Ec 



(12) 



(13) 



In a strictly steady state, dB/dt = identically. How- 
ever, the magnetic flux threading the disk interior to 
a fixed radius r can change on the accretion timescale 
due to advection of field lines by the accretion flow and 
diffusion of fleld lines through the non-ideal disk mate- 
rial. I therefore consider disk solutions that are steady 
on timescales longer than the dynamic time r/VK but 
relax the steady-state assumption on timescales compa- 
rable to the accretion time ^ f / 1 1 ■ 
Expanding Equation (jl3p yields 



dBr 
dt 



dB, 



dz 



B 



B. 



dz 



dVr 

dz 
dEccf) 
dz 



dBr 
dz 



dB^ 
dt 



— Br 



dVs 



dr 



Bz 



B, 



dVr 

dr 



dV^ V^Br 
dz r 
dB. 



-Vz' 



dr 



+ c- 



dEc 



d_B^ 

dz 
dE, 



-B, 



dr 



dz 



(14) 



dV^ 

dz 

(15) 



and 



dB, ^VzBr VrBz 

dt r r 
dBz 



Vz 
dVr 



dBr 
dr 



B, 



E, 



Bz^ C 



dr 



dr 



'C0 

r 



dVz_ 
dr 
dE, 



C0 



dr 



(16) 



where the solenoidal condition (Equation (ITT]) ) was used 
to obtain Equation p^ . Now, ~ Vk, \Vr\ ~ Cg 
through much of the disk and \Vr \ ~ Vz near the top of 
the disk (see Figure 2), z/r < /ix/f = Cs/Vk <C 1 (where 
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Ht is the tidal scale height), and cEc ^ V x B. These 
relations imply that the terms Vr dB^/dz, Bz dVr/dz, 
Vz dBr/dz, Br dVz/dz, and c dE^^jdz in Equation p4l) 
describe processes operating on the dynamic timescale. 
Similarly, the terms B^. dV^/dr, Bz dV^/dz, B^/r, 
Vz dB^jdz, B^ dVz/dz, c dEcz/dr, and c dEcr/dz in 
Equation dTSl) are associated with processes operating on 
the dynamic timescale. However, every term in Equation 
(fT6)) describes a process that operates on the longer ac- 
cretion timescale. I therefore ignore dBr/dt and dB^/dt 
in the r and components of the induction equation, 
since those terms represent changes in the magnetic field 
on the accretion timescale and are negligible compared 
with the terms representing changes in the magnetic field 
on the dynamic timescale, but I retain dB ^Jdt in the z 
comp onent of the induction equation (cf. iKonigl et al.l 

Mm- 

The components of Equation (fT2|) thus become 



dz 



0, 



dEr dEz 
dz 



dr 



and 



dBz 
dt 



■djrE^) 
dr 



(17) 



(18) 



(19) 



The r and z components of the induction equation to- 
gether imply that E^ is constant with height at a given 
r but need not vanish identically. 

2.1.2. Self-Similarity Equations 

The disk equations reduce from partial to ordinary dif- 
ferential equations under the assumption of radial self- 
similarity, which implies that all quantities are described 
by power laws in the spherical radius R for a fixed in- 
clination angle. I normalize all quantities by midplane 
values on a reference magnetic field line [sub script '*'), 
using a scheme similar to that found in O (|1995f ) and 
L96. Magnetic field lines can be labeled by the magnetic 
flux function which satisfies V^* x (p — Bpr, where 
Bp is the poloidal magnetic field. No te tha t the notation 
used here correspo nds to that of iLij (|1995f) and L96 but 
differs from that of IKonigl etaUpOlfj} . where the flux 
function is labeled A and ^ is used to denote the poloidal 
flux; the two quantities are related hy A = 5'/27r. 

Under the self-similarity assumption, the midplane val- 
ues of all physical quantities scale as power laws in ^E*, 
and 

^ cx r^/^^ . (20) 

The magnetic flux distribution parameter C is denoted 
by C in L96. The cold-wind model of BP82 has C = 4/3. 
In a disk that drives a magnetocentrifugal wind, C < 4/3, 
and if the outflow rate in the wind is much smaller than 
the accretion rate, as is the case in protostellar disks, 
then <5C = 4/3-C<C (L96). 

In terms of the self-similarity dimensionless variables 
if = 5'/^'* and s = z/r (denoted by t in L96), various 
physical quantities are normalized as follows: 



r = r^.^p'^ x (s) ; 



2-3C ^ 



P{s) ; 



z = r^ip'^ sx (s) ; 



V = Vk*ip 2 v{s) : 



B = Bz.^'-^^ b{s); 



E 



Vk*Bz* i_ 
V 



e(s): 



Vk*Bz* i_ 



ec (s) ; 



47rr*Vk^ 



-if 



6r(s); 



cB 



Anr 



3 (s) ■ 
(21) 



The variable raised to various powers converts mid- 
plane values on the reference field line (subscript '*'), 
which is labeled by flux function value 'I'*, to midplane 
values on a generic field line (subscript '0') with flux func- 
tion value ^ = ^'q. The combination r^,(p'^ , for example, 
is equal to rg, the midplane cylindrical radius of the field 
line labeled by flux function value V^q. The exponents for 
the power laws in ip in Equation (1211) come from dimen- 
sional analysis of the governing equations. The variable s 
expresses position along a field line. The new coordinates 
{s , if} are expressed in terms of the spatial coordinates 

i/C 



{r , z} as s = z/r and ip = (77^(^77)) 



This gives 



V dr 



dip 
dr 



(r 



+ x's 



ds 
dz 

dip 
97 



Cr 



(22) 



In these relations, x' = (^) 



In terms of the new 



coordinates, the spatial derivatives of a gene ric qu antity 
A = A*(p"''a (s) take the following forms (|Lilll995[ ): 



dA 
dr 



dA 
ds 



dA 

dip 



d(p 
dr 



Af.ip"" (x -I- x's) — C,a'xs 



r^ip"- 



(23) 



and 



dA 

dz 



dA 
ds 



dA 

(^a'x — OiaCLx' 



dip 
dz 



r*ip' 



(24) 



where a' 



^ ds J ip' 



The partial differential equations describing the disk 
can now be recast as ordinary differential equations and 
algebraic relations in the normalized quantities. The 
equation of mass continuity becomes 



, I Vz - VrS ^, / 5 2\ 
■y, - sv^ H z p = 1 2 ~ ) ■^'^ 

2 7\ x' 

^ - 2 ) — ("^ ~"'-^)' (25) 

and the conservation of momentum (Equation ([5])) be- 
comes 

i^z - VrS)v^ ~ —p = y + f0 - (W^ - VrS) — 



1 



(l + s2)3/2 



1 + 

X 



in the r direction, 

in the 6 direction, and 



VrVih , . VAX 

— {Vz - VrS) 



2x 



^ (jzbr - jrbz) 
P 



(27) 



[Vz - VrS) Vz + -P = [Vz - Vr-S) — — 

pi ZX 

2\ x' vx 
- ( 4 - - j 0- + y (>5^ - ubr) (28) 



(1 + .s2)3/2 



in the z direction, where 

^ cl/vl 



and 



_ 2 / 2 



(29) 
(30) 



The quantity v\ q/v^ is denoted by a in L96. 
By definition, \I> = ^^ip. The relation = - (||)^ 

therefore gives {Bz^ip^^^'^br) {r^ip'^x) = ^^(^)- 
From the definition of the flux function, ^ = Bz*, so 
we have 



c' = br 



(31) 



At a point (r, z) = r*(p^ {x{s), sx{s)) on a field line, the 
unit tangent vector along the field line is 

x'f + {x + x' s) z brf + bzZ 



+ {X + X'sf VW+^ 
This relation, combined with Equation pip , gives 

bz = brS + \ . 

x'^ 



(32) 



Equation (1321) . which expresses bz in terms of b^, x, and 
s, automatically satisfies the solenoidal condition (Equa- 
tion (fTTI) '). 
Ampere's law (Equation (|9])) becomes 



bL 



in the r direction, 



1 \ 

- - 2 1 —60 - xjr 



b'^ + sb'^ = xjif, + \^-~2 
in the 6 direction, and 



— {br + bzs) -I- bz 

X 



(33) 



(34) 



(35) 



in the z direction. Using the algebraic relation for 6^, the 
(26) (j) component simplifies further to 



b' = 



+ (i-2) ^^{br+bzS)+b 



1 + 



(36) 



The poloidal components of Ampere's law (Equations 
(|33)) and (|35|) ) together imply 



J.=>s+Q-l)^. (37) 
Ohm's law (Equation ^) becomes 



J = cr • Cc = (Toecii + ^nb X ec_L + ffpeci , (38) 
where the comoving electric field is given by 

ec = e + v xb (39) 
The three components of Equation ([55)1 are 

j,. = aoScrW + -y {b^ecz± - bzCc^^) + apecr± , (40) 



j<p = ^oec<^|| + — {bzecr± - brecz±) + ^pec(/.± , (41) 



jz CToeczll + — {brec4,± ~ b^ecri_) + ^peczi_ , (42) 



where b— \J b'^ + b'^ + bl , and the magnetic field-parallel 

and -perpendicular components of the comoving electric 
field satisfy 



and 



br b^ bz 
ecr±br + ec0±60 + ecz±bz = , 

ec = ecii -I- ec± . 



(43) 
(44) 

(45) 
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The poloidal components of the induction equation 
(Equations pT)) and p9|) ) yield the relation e^{r,z) = 
60 (r, 0), which implies 



(46) 



The (f) component of the induction equation (Equation 
pS])) yields 



(47) 



Combining equations (|25]) . p6|. and p8| yields 



(w^ - ^;^s)^ - 61 (1 + s2) 



2- - 



(Wz - VrS) 



+ 



(w2 - t;rs) 



+ 4 



(j0 i^- 



vis 



bcf, [ir+jzS]) 



(48) 



This relation shows that the point where [v^ — v^s) = 
(l + s^) is a singular point of the system of equa- 
tions, since the denominator of p' vanishes there. The 

quantity \J {v^ — Vj-sf' / (1 + s^) is the component of the 

poloidal velocity perpendicular to the local surface of self- 
similarity (the surface of constant s), normalized by the 
Keplerian velocity. The denominator of p' vanishes when 

{vz — Vrs)^ / (l -I- s^) ~ 9 — (cs/Vk)^. This occurs when 

the unnormalizcd velocity component Vz — Vrs/ \\ + s^ 
becomes equal to the local sound speed; this defines 
the sonic point. The Alfven point is similarly a critical 
point of the cold ideal wind equations, and occurs where 
the same poloidal velocity component Vz — Vrs/^/l + 
reaches the local Alfven speed. 

Combining Equation p8]) with Equations ((26| and ([28]) 
gives expressions for v'^ and v'z- Equation pip gives x' , 
Equation (f27|) gives u^, and Equations (|36l) and (pS)) give 
6J, and 6^, respectively. The final ordinary differential 
equation is Equation (|T7)) . Solving these 8 ordinary dif- 
ferential equations requires knowledge of the conductiv- 
ity profile (To(s), (Th(s), and a'p(s), which can be ob- 
tained using a known disk ionoization profile or some 
simplifying assumptions. In order to solve for ej, in Equa- 
tion (|T7)) . it is necessary to express Cz as a function of 
the other variables and then differentiate the resulting 
relation. This expression is derived in Appendix \K\ 

2.1.3. Values at the Midplane 

All physical quantities are assumed to have either odd 
or even symmetry in z about the midplane. The mag- 
netic field is vertical at the midplane, so the s derivative 
coincides with a z derivative at that location. 

By definition, xq = I and po = 1- Since the magnetic 
field is vertical at the midplane, &r,o = and b^^Q = 
0, and since the vertical component of the flow velocity 
must be continuous across the midplane, Vz o = 0- The 



midplane radial velocity and azimuthal electric field are 
given by 



and 



EB = -cE^fl/CsBz 



(49) 



(50) 



In dimensionless form, these relations are Wj.,o = —(-B^l"^ 
and 60.0 — — eB^*^/^, which combines with Equation (|46l) 
to give 



i _ 5 
i 2 



(51) 



The quantities e and eb arc parameters of the model 
whose values are set by conditions that must be satisfied 
at special points in the flow (see Section [Sj. The deflni- 
tions of e and ee given here are the same as in WK93; a 
factor of 0^/2 appears in the expressions for w^.o and e^^o 
because velocities here are normalized by wk rather than 
by Cs as in WK93. 
At the midplane, Equations and (|27l) become 



2 



0,0 







and 



+ Vjrfi = , 



(52) 



(53) 



respectively. Since the magnetic fleld is vertical at the 
midplane, 6cr-||,o = 0, ecr±,o = Ccr.o = e^.o + W0,o, 
^0011,0 = and 6c0j.,o = ec0,o — e<p,o — Vr.o — 
{e — es) 0^^^ . The radial and azimuthal midplane cur- 
rent densities are therefore jr.o = — o'H,oec0,o + o'p a^cr.o 
and jcf,.o — o'H,oecr,o + ^p.,oSc<j>,o (see Equations ([to|l and 
((4T|) ). Substituting these relations into Equations (|52p 
and ([55)1 gives 



1 



( CTH^o [e,.,o + W0,o] + CTp^o [e - cb] ) = , (54) 



and 



e0i 



W0,O 



-o-H,o e - £B 



+ CTp,o [er,o + V0,o]) = . (55) 



Combining these two equations in the unknowns v^^ and 
Gr,o gives 



2 , / ^H.pe^^ 



- 1 



4- 



'H,0 ~r "P,0 
O'P.O 



(e-eB)0^ 



0, (56) 



which can be solved for w^^o- Either Equation ((54)) or 
Equation (|55p can then be used to calculate e^o- 
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TABLE 1 

Summary of Midplane Values 



Variable Midplane Value s-Derivative At Midplane 

X 1 

p 1 

Vr —€0 2 

Vz 



e02 



5 _ 2 

b<4 -ir,0 

bz 1 



-eB6'2 



The values of several physical quantities and their s 
derivatives at the midplane as functions of the model 
parameters and i^^^o s^re summarized in Table [TJ 

Note that the expression for v'^ q is calculated using 
the mass continuity Equation ()25p . rather than Equation 
which is identically at the midplane. Also, from 
Equation (PSjl we have jz.o — 0. The magnetic field is 
vertical at the midplane, so Ccz.o = Cc^n^o- From Ohm's 
law (|42)) we have jz_o = <5'o,oecz,o; since ao,o is non-zero, 
jz,o = implies ecz,o = 0. 

The midplane is a singular point of the differential 
equations describing the disk, and special care is needed 
in initiating the numerical integration from that location. 
A common choice for dealing with this issue is to approx- 
imate the values of the dependent variables a small dis- 
tance above the midplane using Taylor series expansions, 
a task that is made much easier by the assumption that 
all physical quantities have either even or odd symmetry 
about the midplane, which implies that the Taylor series 
for all quantities contain only even power terms or only 
odd power terms — for instance, the fact that 6J. g ^ 
immediately implies that 6"o = 0. Looking at Table [U 
it is clear that second-order Taylor series expansions for 
Vz, br and 6^ only have one term apiece, corresponding 
to the first derivatives of those quantities at the mid- 
plane, which can be calculated using the expressions in 
Table [T] and Equation ([55]) for w^ q- The second-order 
Taylor series expansions for x, p, Vr, v^, and include 
terms involving Xq, pg, ^"g, w^'g, and e^g; the procedure 
for c alculating these quantities is described in Appendix 

EH 

2.2. Amhipolar Diffusion Case 

The ambipolar conductivity regime occurs when all 
charged species j have j3j ^ 1. Substituting into Equa- 
tions jSj-lISI) gives (To ^ o'p ^ I^hI in the ambipolar 
diffusion limit. 

I assume that the charged species consist of ions (sub- 
script 'i') and electrons (subscript 'e'). Charge neutrality 
and a comparison of the rate coefficients for electron- 
neutral and ion-neutral momentum exchange given by 
iDraine et al.l ([T993i) imply that the heavier ions dominate 



the charge-neutral drag force. Combining this with the 
balancing of the charge-neutral drag and Lorentz forces 
yields 

^^=7Pi(^i-V) , (57) 
cp 

where is the ion velocity and pi is the ion mass density. 

For a disk in the ambipolar diffusion limit where the 
only charged species are ions and electrons, the Elsasser 
number is A = ^pi/Q.]<^ = T, the ratio of the dy- 
namic timBscaIe_to_thfijier^^ momentum exchange 
time (|Konigl fc Salmeronll201 1'). This ratio is denoted 
by T] in L96. The neutral-field coupling parameter T 
in this case is related to the Pedersen conductivity by 
T = {4Trvl/c^){ap/nK) (SKW07). 

The tight coupling of the ions to the magnetic field 
in the ambipolar diffusion limit implies that the comov- 
ing electric field vanishes in the frame of the ions. The 
electric field in the disk is therefore 



E = --Vi X B. 



(58) 



2.2.1. Governing Equations 
Equations and (1551) yield the induction equation 

dB 

'dt 



V X (Fi X B) 



(59) 



From Equation (|58|) and the definition of ee, cE^ — 
VirBz - VizBr = -eBCsir)Bz{r,z^ 0), or 



VizBr = VizBz + eBCsBz(r, z = 0) 



(60) 



where all quantities [except the midplane quantity 
Bz{r,z = 0)] are evaluated at a point (r, z) on the field 
line of interest. 
The (p component of the induction Equation (|12p is 



\ dr dz ' 



(61) 



or 



^ (y.^Br - V,rB^) = — (VizB^ - %Bz) . (62) 
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2.2.2. Self- Similarity Equations 

The ordinary differential equations obtained under the 
radial self-similarity assumption are identical to Equa- 
tions (12)-(14) and (16)-(21) in L96. The equations 
used here differ from those in L96 only in the treat- 
ment of the induction equation; I allow for a nonzero 
ee (or, equivalently, a nonzero Erf,)^ which changes the 
usual flux-freezing condition for the ions. As a result, 
Equation (10) in L96 is replaced by Equations (15(1)) and 
((^ . Equations (22)-(24) in L96 (the self-similar, nondi- 
mensionalized versions of Equation (10) in L96) are re- 
placed by the self-similar, nondimensionalized versions of 
Equation (pU)) . 



and Equation ([S^ . 



(63) 



brX^ {Virb^, - Vi^br) , (64) 



where the normalization scheme from Section 12.1.21 is 
used, with the addition of 

Vi = WK*<y5"^l'i (s) ■ 

Taking the dot product of Swith both sides of Equa- 
tion (157)) gives B- iVi — V) — 0, which results in the 
following algebraic relation for Wj^: 



V{^b^ ^brVr + b^v^ + b^v^ - V-ir [br + y 

_eB(?l/2^^1/C-5/2^ 
br 



(65) 



which can be differentiated to give v[a for use in Equation 



I now have a system of 8 ordinary differential equa- 
tions in the variables x, p, Vr, v^, Vz, fir, br and b^. This 
system differs from that of L96 in the treatment of the 
induction equation. A solution of this version of the disk 
equations and a matching wind solution are presented in 
Section |4l That solution assumes a functional form for 
the matter-field coupling parameter T oc pi that differs 
from the one used in L96 — while L96 assumes that T 
increases monotonically with height, I take T to be con- 
stant along a field line, as would be the case when the 
ion density is constant (WK93). In fact, the coupling 
parameter is expected to decrease with height in the up- 
per reaches of real protostcllar disks not too far from the 
origin (iSalmc ron &: Wardle, 2005.) , reflecting the decrease 
of the ion mass density pi with height (which occurs even 
though the ionization fraction increases with height be- 
cause the neutral mass density p drops rapidly with z). 

On the other hand, the coupling parameter increases 
with radial position because the Keplerian frequency JIk 
decreases. Taking T constant is a compromise between 
these competing effects as the integration proceeds up- 
wards and outwards along a field line. Finally, the solu- 
tion presented here differs from the solutions presented in 
L96 in that I initiate the numerical integration of the disk 



equations using a Taylor series expansion (see Appendix 
IB. 2)) . whereas L96 uses the approximation t;^ = u^s in a 
small region near the midplane to obtain algebraic rela- 
tions for Vr and v^. 

2.2.3. Values at the Midplane 

From Equation ()63l) . v\rfi — —£bS^^^- The r com- 
ponent of momentum conservation (Equation ([8])), with 
Equation ()57p used to replace the Lorentz force term with 
an ion-neutral drag force term, yields Equation (13) of 
L96. At the midplane, this equation reduces to 



2 



-1 



4- 



' + To {Virfl ~ Vrfi) = , (66) 



or 



<o = l- — - (^4--j(? + To^^(e-eB) , (67) 

which gives Vrp^ in terms of the parameters 9, ^, Tq, 
e, and e^. Continuity of fj^ implies vi^fi = 0, and the 
(f) component of the momentum equation (cf. Equation 
(14) of L96) evaluated at the midplane gives 



To {Vict^fi - Vrjj^o) 



(68) 



which can be used to calculate Vi^fl. 

Combining Ampere's law (Equation ([9|)) with Equation 
(|57)) yields equations for b'^ and b'^ (cf. Equations (17)- 
(19) of L96); evaluating at the midplane yields 



and 



00,0 = (^^i0,O - V4,.o) 



(69) 



(70) 



The expression for v'^ q (see Tabled]) remains unchanged. 

The Taylor expansion procedure used to initiate the 
integration is described in Appendix IB. 21 

3. MODEL PARAMETER CONSTRAINTS 
3.1. Model Parameters 

The dimensionless parameters 9 = c^/v^ and ly = 
v\ o/^K model inputs, as are the conductivity tensor 
component profiles ctq (s), an (s), and ap (s). Alterna- 
tively, values for 6 and together with a profile T (s) for 
the neutral-field coupling parameter suffices for the case 
of a disk in the ambipolar diffusion regime. In either case, 
a disk/ wind solution involves constraining three param- 
eters to satisfy certain conditions imposed by physical 
processes outside the disk. The first constrained param- 
eter is C, the power-law index of the magnetic flux distri- 
bution (Equation ()20p ). The second is e, the normalized 
midplane value of Vr (Equation (HHj))- The third is ee, 
the normalized midplane value of (Equation ()50p ). 

Two of the conditions to be satisfied are that the 
disk solution pass smoothly through the sonic point and 
match onto a wind solution that passes smoothly through 
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the Alfven point. The third condition relates the angle of 
the magnetic field at the disk surface to the distribution 
of the magnetic field exterior to the disk (see Section l373)) . 

3.2. Matching to a Wind Solution 

The sonic point (z = Zgp) is a singular point of the sys- 
tem of equations describing the disk (see Equation (|48l) 
and subsequent discussion), so it is not possible to inte- 
grate directly into the sonic point, but it is possible to 
get close. I run a bisection routine on e, holding ^ and 
eb fixed, until I obtain a disk solution that accelerates 
smoothly almost all the way up to the sonic point. I then 
extrapolate the solution to Zgp and attempt to match the 
disk solution to a BP82-type wind solution, using the 
procedure outlined in SKWll. That is, I calculate the 
value of the wind parameter k (the normalized mass-to- 
magnetic flux ratio) using the disk solution at the sonic 
point, and I calculate Br^s/Bz,s and As, respectively the 
ratio of the radial and vertical components of the mag- 
netic field (^{, in the notation of SKWll) and the total 
specific angular momentum, using the disk solution at 
the disk surface. Here, the disk surface is taken to be the 
point where the azimuthal neutral gas velocity matches 
the local Keplerian rotation rate. I run a bisection rou- 
tine on ( until the values of k, A and calculated from 
the disk solution match the values at the base of a wind 
solution. 

The matching procedure used here guarantees continu- 
ity of the magnetic field components, the poloidal mass 
flux, and the total specific angular momentum. How- 
ever, it should be pointed out that the electric field is, in 
general, not continuous at the transition from the disk 
solution to the wind solution: in particular, there is a 
non-zero azimuthal comoving electric field Ecp in the disk 
solution because the parameter ee is generally non-zero, 
but the cold, ideal MHD equations describing the wind 
solution have Ec^, = identically. In addition, the BP82 
wind solution is derived for a razor-thin disk that has 
vanishing radial and vertical gas velocities at the base of 
the wind (negligible compared with the azimuthal veloc- 
ity, which is exactly Keplerian at the base of the wind); 
these are good approximations when the region of inter- 
est is the large-scale outfiow, but break down on small 
scales comparable to the thickness of the disk. The wind 
solution is therefore expected to give physically uninter- 
esting results very close to the disk surface. In fact, the 
wind solution has a singular point at the midplane just 
as the disk solution does. I follow BP82 in integrating 
the wind equations downwards starting from the Alfven 
point and stopping very close to the disk surface, where 
the solution is matched to a Taylor expansion. I then in- 
terpolate between the disk solution at the top of the disk 
and the wind solution at a distance x/xa = 0.002 above 
the disk surface, where x = z/rg^ Vs is the radial distance 
of the disk surface along the field line of interest, and XA 
is the value of x at the Alfven point. 

The wind solution extends past the Alfven point. 
There is a second critical point of the equations for a self- 
similar wind corresponding to the modified fast magne- 
tosonic point (BP82). However, before a cold BP82-type 
wind solution reaches the modified fast magnetosonic 
point, the flow "over-coUimates" as the field lines bend 
back towards the axis of symmetry. This unphysical 
behavior results from ignoring heating of the jet ma- 



terial and thermal pressure support in the cold wind 
model. "Warm" self-similar jet models whose solutions 
do cross the modifi e d fast mag netosonic have bee n found 
by iVlahakis et all (pOOOh and iFerreira fc Cassi (2004). 
For this study, I terminate the wind solution when the 
field lines begin to bend back towards the jet axis. 

3.3. The Disk Surface Field Constraint 

OLOl pointed out that in a self-consistent, global 
model of a magnetized disk, the inclination of the 
poloidal field at the disk surface can be determined un- 
der the approximation that the external magnetic field 
is nearly potential (V x B w 0). An analytic expression 
relating i?j.,s to the distribution of B^^s along the disk sur- 
face was derived, under th is approximation and as suming 
an infinitely thin disk, by iCiolek fc MouschoviasI (fT993i) , 
who modeled magnetically supported, se l f-gravi tating 
molecular cloud cores, and bv lLubow et al.l ()1994aD . who 
modeled fiux advection in viscous accretion disks. The 
poloidal field inclination at the disk's surface is directly 
related to the balance between the inward advection of 
magnetic fiux by the accretion flow and the outward dif- 
fusion of magnetic flux through the non-ideal disk. This 
connection follows from the fact that the radial compo- 
nent of the magnetic force (per unit area) at the disk 
surface, which is cx i3r,sSz,s, is balanced in the ambipo- 
lar diffusion regime by the radial ion-neutral drag force, 
which is cx Vir — Vr (or, taking the midplane values, 
cx e — eb, which indicates the explicit relation of 6r,s 
to the disk model parameters). Determining 5^,8 self- 
consistently in this manner generally implies a nonzero 
value of eb- The relevant expression for Br.s is 

rBr,s{r) 

poo pOO 

/ dkkMkr) dr'r'[B,,o{r')~B,,f]Mkr'), {71) 

Jo JQ 

where Jq and Ji are Bessel functions of the first kind 
of order and 1, respectively, B^flir) is the midplane 
vertical field component at radius r, which is assumed 
to remain constant up to the surface [so that Bz,s{r) = 
Bz,o{r)], and By^i is the external magnetic field at "infin- 
ity," which is taken to be ^ Bz^Q{r), so it may be safely 
ignored in the second integral. 

For a power-law profile B^o oc r^". Equation 
dZIl) yields Br,sir)/Bz4r) = [T{1 - a/2)r(l/2 + 
a/2)][r(a/2)r(3/2 - a/2)] for Re{a} G (-3/2,-1/2), 
where F is the gamma function. The power-law exponent 
is related to the disk model parameter ^ by a = 2 — 1/C, 
so 

BrAr) ^ F[1/(2C)] F [3/2 -1/(20] 

BzAr) F[1-1/(2C)]F [1/2 + 1/(20] ' ^ ' 

This result holds for ( S (2/3, 2) and is therefore appro- 
priate for the values C ^ 4/3 that are expected to char- 
acterize physically viable solutions (e.g., L96; F97). For 
C = 4/3, Equation 1^ yields Br^r) / Bz,s{r) = 1.428. 

Once a disk solution that passes smoothly through the 
sonic point has been matched to a wind solution, I bisect 
on EB until the value of Br.s matches the value imposed 
by the field distribution. 

4. ILLUSTRATIVE SOLUTION 
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Fig. 1. — Density structure of the constant T disk solution. 
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Fig. 2. — Neutral velocity structure of the constant T disk so- 
lution. The Keplerian velocity has been subtracted from the az- 
imuthal component, and the plotted velocities are normalized by 
the local sound speed. 

Protostellar disks have typical inferred temperatures 
at 1 AU of - 100 - 300 K, masses of - 0.005 - 0.2Mq, 
column densities at 5 AU of ~ 0.3 — 300 g cm~^, and ex - 
tend to radii - 200 AU fe.g. . I Andrews k WilliaiSl2007D . 
Disk/wind model solutions typically have va,o/cs ^ 1 
(e.g., WK93), and observations of protostellar systems 
with jets indicate that the ratio of the outflow and ac - 
cretion rates is ~ 0.01 — 0.1 (e.g., IHartigan et al.l [19951 ) . 
Finally, observations of remanent magnetization of mete- 
orite samples suggest a magnetic field o f ^ 1 G at 3 AU in 
the early solar nebula (iLevy & Soncttl [l978D . Motivated 
by these facts, I choose parameter values = 1x10"^ 
and v = 7.744 x 10"''. These values together with the 
constant coupling profile T (s) = 10 yield a solution of 
the disk equations described in Section [2^ This solution 
passes smoothly through the sonic point, matches to a 




1 ^ I 

0.05 0.1 0.15 

Fig . 3 . — Magnetic field structure of the constant T disk solution. 
All components are normalized by the midplane field strength. 
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Fig. 4. — Poloidal magnetic field line shape of the constant T 
disk solution. The location of the sonic point is marked ('*')■ 

wind solution that passes smoothly through the Alfven 
point, and satisfies the constraint on the magnetic field at 
the surface of the disk. The corresponding disk and wind 
parameters are S( = 3.05 x 10~^, e = 0.45, ee = 0.03, 
K = 1.27 X 10-3, A = 108.4, and = 1.42. The small 
value of SC found here is consistent with the predictions 
of L96 and F97 that cold wind-driving disks should have 
field-line distribut ions very close to the h miting case of 
C = 4/3 (but see iCasse fc Ferreiral [2000l for arguments 
that larger values of SC can occur in wind-driving disks 
when thermal effects are consideretfl) . 

^ T he quantity is related to the quantity § in lCasse & Ferreiral 
l |2000 ) as (5(^ PS 8/9 J; the ratio / of the mass outflow in the wind 
between radii rjn and rmax to the mass accretion rate at rin is 
/ = (rmax/nn)* " 1- For (rmax/nn) ~ 1000 and an observed 
/ 0.01 - 0.1, 5 is inferred to lie in the range g ~ 0.0014 -0 -014. 
Of the two representative solutions in lCasse &: Ferreiral l|2000l ). the 
one with almost no entropy generation gives ^ ~ 0.001, while the 
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Fig. 5. — Density structure of the constant T disk solution and 
the matched wind solution. 
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Fig. 6. — Velocity structure of the constant T disk solution and 
the matched wind solution. The plotted velocities are normalized 
by the local Keplerian velocity. 

A rough estimate of the ratio of the mass outflow rate 
from both sides of the disk, between an outermost radius 
f'max and a given radius r, to the mass accretion rate 
through the disk at r, can be derived from Equation (89) 
of SKWU as Mwind/Mdisk « ln(w/r)/[2(A - 1)]; for 
the derived solution, and taking rinax/?" ~ 1000, this ratio 
is « 3.2 X 10~^. For a central protostar of mass 0.5 Mq, 
the model has a column density of « 42 g cm~^ at 1 
AU and a column density of « 19 g cm~^ at 5 AU if 
the density is normalized so that the accretion rate is 
« 3 X 10^^ Mq yr~^ and the total disk mass out to 200 
AU is « 0.055 M0. In compariso n, the minimum mass 
solar nebula model ()Havashilll981f) has a column density 

one with strong entropy generation gives § ~ 0.5, in which case 
almost all of the disk mass flows out into the wind rather than 
accreting onto the central protostar. 



Fig. 7. — Magnetic field structure of the constant T disk solution 
and the matched wind solution. All components are normalized by 
the midplane field strength. The radial and toroidal components 
are at the midplane, and the vertical component is 1 at the 
midplane, but all three components change rapidly within the disk, 
which extends up to z/ro ~ 0.079. 
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Fig. 8. — Poloidal magnetic field line shape of the constant T 
disk solution and the matched wind solution. The location of the 
Alfven point is marked ('*'). 

at 1 AU of w 1700 g cm-2 and at 5 AU of w 150 g cm^^. 
The density normalization used here corresponds to a 
midplane magnetic field strength at 1 AU of |B| « 2.5 G 
and at 3 AU of « 0.6 G. The disk temperature at 1 AU 
is « 125 K. 

The gas density, velocity, and magnetic field struc- 
tures of the disk solution as functions of position along 
a field line are shown in Figures [TH31 The plots ex- 
tend from the midplane, z = 0, to the sonic point, 
2sp ~ 0.14 ro (~ 3.8 /it)- The location of the sonic point 
is marked in Figure 01 which shows the poloidal magnetic 
field line structure in the disk. The magnetic field com- 
ponents and the density are normalized by the midplane 
values of B and p, respectively. The velocity compo- 
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nents, with the Keplerian velocity subtracted from the 
azimuthal neutral and ion velocities, are normalized by 
the local value of the sound speed along the field fine 
(that is, taking into account the spatial variation of the 
sound speed). 

At the midplane, the field is vertical by symmetry; 
above the midplane, the radial component grows rapidly, 
a small negative toroidal component develops, and the 
vertical component decreases slightly. The vertical veloc- 
ity Vz starts at zero at the midplane, decreases to small 
negative values in the near-midplane region (refiecting 
the radial convergence of the accretion flow), and then 
increases back through zero at 2; « 0.054ro(~ 1.5 /it), 
becoming comparable to Cg at z w 0.112 tq (~ 3.1 /it)- 
The radial velocity starts out at — 0.45cs at the mid- 
plane and decreases in the near-midplane region, reach- 
ing a minimum at z 0.016 tq (~ 0.5 /it) and then 
increasing, switching from an infiow to an outfiow at 
z « 0.071 ro(« 2 /it). The azimuthal velocity is sub- 
Keplerian at the midplane 0.93 Vk.o) but increases 
smoothly with height, becoming super-Keplerian at the 
disk surface, Zg ~ 0.77ro(« 2.2 /it), slightly above the 
point where Vr becomes positive. In the hydrostatic 
approximation {Vz =0), Vr = and = Vk at the 
same point (WK93). The effective density scale height 
/loff (defined so that p (z = hcs) / p (z = 0) — l/-\/e) is 
hcff w 0.22 /iT, indicating that magnetic squeezing has a 
significant effect on the disk structure (WK93). 

The gas density, velocity, and magnetic field structures 
of the matched disk and wind solutions as functions of 
position along a field line are shown in Figures [5H71 The 
plots extend from z = to z « 48 tq, where the field lines 
begin to bend back towards the jet axis and the wind 
solution terminates. Along the way, the wind solution 
passes through the Alfven point, located at z « 38.9 tq. 
Figures [5]-[7] are shown with the variable z/tq plotted on 
a logarithmic axis so that the structures of the disk and 
wind regions can be distinguished. The location of the 
Alfven point is marked in Figure [8l which shows the 
poloidal magnetic field line structure in the disk and 
wind. In the region above the disk, the poloidal field 
quickly drops while the azimuthal field amplitude first 
grows slightly (S^ becomes more negative) and then de- 
creases, in accordance with the expected behavior of a 
wind-driving disk. The vertical and azimuthal velocities 
increase smoothly throughout the outfiow, while the ra- 
dial velocity reaches a maximum and begins decreasing 
well below the Alfven point. Note that Figure [5] differs 
from Figure [5] in that the unmodified azimuthal velocity 
is plotted and the velocities are normalized by the local 
value of the Keplerian speed (that is, taking into account 
the spatial variation of the Keplerian speed). 

The solution presented here is similar to the WK93 
model in that it considers a strongly coupled, ambipo- 
lar diffusive disk, but differs in that it is a global model 
with a non-zero value of eb determined self-consistently 
from a constraint on the field line inclination at the disk 
surface. The equilibrium solution curves for the present 
model are undoubtedly modified from those of K04, but 
on the assumption that the changes are not drastic, a 
comparison of the representative solution and the results 
given in K04 may be instructive. The representative so- 
lution given here has k = 1.27 x 10"'^, A = 108.4, and 



1/2 

a = {(t/9) ' = 0.88; the quantity ry in K04 corresponds 
to T. Comparing these values to Figures 1 and 2 of K04, 
it appears that the representative solution given here lies 
near the lowest point of the lower branch of the rj ~ 10 
equilibrium solution curve in Figure 1, and to the left 
of the upper branch of the rj = 10 equilibrium solution 
curve in Figure 2. In both cases, the comparison suggests 
that the representative solution given here is stable. 

This illustrative solution assumes that the underly- 
ing disk is in the ambipolar conductivity regime. Am- 
bipolar diffusion is likely a reasonable assumption in the 
outer regions of typical protostellar disks, but may not 
hold throughout the inner regions of protostellar disks, 
which are more likely in the Hall or Ohm conductivity 
regimes. However, it is worth noting in this connection 
that models of wind-driving disks have column densities 
in the inner regions that are much smaller than in typical 
turbulent-viscosity accretion disks. These smaller col- 
umn densities suggest that the region of a wind-driving 
disk in which ambipolar diffusivity is a good approxima- 
tion may be substantially larger than in other types of 
disk models. 

5. CONCLUSION 

The model for the global structure of a magnetized pro- 
tostellar disk/ wind system presented here improves on 
previous studies (e.g., WK93; L96; F97; OLOl) by includ- 
ing, for the first time, a fully self-consistent treatment of 
flux migration on the accretion timescale. This treat- 
ment employs the approximation that the magnetic field 
on scales z < r above the disk surface is very nearly po- 
tential, leading to a simple constraint on the surface field 
inclination (Equation (ITT]) '). The model also improves on 
previous studies (e.g., SKW07) by extending the con- 
ductivity tensor formalism for describing non-ideal disk 
material from radially localized disk models to global, 
radially self-similar disk models. 

The model presented here describes solutions for a 
disk and wind that are steady on timescales up to the 
flux migration timescale. Studi es of simplified equ i- 
librium mode l s of d isk winds by 'Lubow ^ al.l (jl994bfl , 
ICao k Spruit (|2002D . and [Campbell (200"9f ) have raised 
the possibility that magnetocentr ifugal disk winds are in- 
herently unstable. iKonigll ()2004l . hereafter K04) showed 
that for strongly coupled wind-driving disk models of 
the sort described in WK93 there are two equilibrium 
solution branches, one stable and one unstable, and ar- 
gued that protostellar systems correspond to the stable 
branch. This line of argument has received some support 
from simulations of axisymmetric, diffusive accretion 
disks launching self-collimating, magnetocentrifugal disk 
winds that have been carried out by Cassc & Kcppcni 
(2002", '2004) and lZanni et al. I ((200711 . These simulations 
confirm that the resulting large-scale jets are continu- 
ous and steady over many dynamic timescales, provid- 
ing support for the present model's quasi-steady state 
assumption. In turn, the semianalytic model presented 
here facilitates calculations that cannot currently be car- 
ried out in full numerical simulations. 

The formulation presented here allows construction 
of global disk/wind models with realistic conductivity 
profiles in the disk. Finding matched disk/wind solu- 
tions of the sort presented here is a difficult task due 
to the imposed restrictions of requiring all angular mo- 
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mentum transport to occur via the wind, and the strict 
self-similarity of the disk and wind solutions. In real- 
ity, protostellar disk/wind systems may transport angu- 
lar momentum simultaneously through multiple mecha- 
nisms, such as magnetic braking, gravitational instabil- 
ity, and turbulence (see SKW07 for a radially localized 
disk model that includes angular momentum transport 
by both a wind and turbulence driven by the magnetoro- 
tational instability), and are unlikely to be globally self- 
similar. However, simple models of the sort presented 
here are valuable tools for investigating global processes 
in disk/wind systems, such as dust transport and pro- 
cessing (S. Teitler, in preparation). Studies of dust evo- 
lution processes in wind-driving disks will complement 
similar stu dies carried o ut in simple models of turbulent 
disks (e.g., iKeller fc Ga il 2004; Brauer et al. 2008) and 
in the X-wind model (e.g., iShu et al.l 120011) . and are of 
particular importance for wind-driving disks since small 
dust grains can have strong effects on the coupling of 



neutral disk materi al to magnetic fields fe.g.. lNishi et al.l 
flMtlWa^da [20071) . 
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APPENDIX 

CALCULATING E'^^ 

Given the values of the independent variable s and the dependent variables x, Vr, v^, Vz, br, b^, e^, and p, the 
auxiliary variables b^ and can be obtained directly from the algebraic relations (|32|) and (|46|) . Their derivatives are 



b' = b's + br 



2^ =fe;s-6, 
x-^ 



and 



60 (0) brX< 



The auxiliary variables Ccr and ec0 are obtained from Equation (|39|) . 
The next step is to combine the algebraic relations (|57)) . PO]) . and (H^ to obtain 



o'oeczii + {brec^± — b^ecr±) + <^pecz± = ^^oecrli + (b<}>^cz± — ^zGc^x) + ffpecr l I s 



Using Equations and (03]), we have 



^cr 6cr±; 



1 



^ - 1 



C J X 



(Al) 



^c^ll — 7 \^Q,r 6cr_Lj 7 ^c0_L — , l^^cr ^cr_L j ? 



_ ^ / _ \ 
Cc0 ^ \^Q.r Ccr_L J 



(A2) 



giving e(,|| and ec_L in terms of e^rL and the known quantities Ccr, &c<^>^ and h. Substituting these expressions into 
Equation (jAip and solving, we obtain 



6cr 



br I brbz b 



bri, + S 



+ 



bd,b. 



brbz b. 



•(POz 




b. 



+ br 



i_Ab_, 
C J X 
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which simphfies to 



where 



and 



B 



ecr± = ACcr + B = ACr + B + A {Vrb^ — V(f,hr) , 

dobbz + a-px^bb'^ — sanx'^b'^b^ — (Jabrb, 

A = 



[('''O ^ '''p) bbz + apx'^b^ — sa-nx'^b'^b^] 
{—dpx'^bbrbffi + sanx'^b'^br + (TH^r) ^ ^ xbbrb^bz 



[((To — <Tp) bbz + apx^b-^ — saux^b'^b^ 



(A3) 



Equation (jA3p for ecr± in terms of known quantities can be used to infer ecr||, ec0||, ec0±, ecz\\, ecz±, Scz, Sz, jr, j^, and 
jz using Equations (jA2l) and (|40)) -(|42 |) . These auxihary quantities can in turn be used in Equations (l48l) . ([26)) -(l28 | . 
(|36)) and (|33l) to obtain p', w^, u^, w^, 6^, and fc^. 
The differential Equation (|47p for involves the derivative of the auxiliary quantity e^; is given by 



= Sczll + GczX + Vfi,br - Vrb^ 



cr *=-cr_ 



^^c0 ^ (Ccr ^^cr_L) 



+ V^br — Vrt 



Upon substituting for ecr± from Equation (|A3|) . this becomes 



\br bz brbz J 

+(..5z-.zMf^i-A]-^A+l-[i-A; 



-B- 



brbz 



Differentiating Equation (IA4[) for Cz gives 



where 



-ec,/, + Vrpbr - Vrbff, ■ 



e'z = Ce'^^D, 



(A4) 



(A5) 



C=^(l-A)-^^+-f (1-A) = ^ + 



bl 



D = e'er + C (v'^bz + V4,b'z - v'zb4, - Vzb'^) + C {v4,bz - Vzb4,) 



-B 



brbz 



2B- 



bb' 

brbz 



B 



. bX 

b^bz 



B 



ybi 



brbl bz 



6? 



-^^0 + + v^b'^ - v[.b^ - Vrb'^) 



brbz blbz brbl 



blbz 



brbl 



^ \o'Qhhz + (Tob'bz + a-obb'z + (TpX^bb^ + 2apxx'bb^ + apx^b'b^ + 2apx^bb^b'^ 

2 1,2 1 ~l 2i 

—crux b 0^ ~ S(Tp^x t 



- ^ [aobb^ + a-px^bb^^ - 
+ {do-dp)bb' ' — -^'^3 



sa-'j^x'^b^bcf, — 2sd-iixx'b^b^ — 
- ^nbrb'^] ^ y/ [((5-0 - 

4 - s&ux'^b'^b^ - aiibrb^] [ (ctq - 5-p 



^- . . 'b: 

'd'pibrb^ ~ anb'rb^, ~ dYibrb'^] ^ j [((Tq - o-p) 662 + apa:^6^ - so 



2saHX^bb'b^ ~ sanx^b^h 

saiiX^b^bif 



Tp) bb^ + (cto - CTp) 
o-px^fe^ + 2apxx'b^ + Sapx'^b'^b' - a^x'^b'^b^ ~ sa'p^xHH^ 
2s(JHXx'b^b^ — 2sd]iX^bb'b(j, — sd^iX^b^b'^i^ j [(cto ~ ctp) bb^ + f7px^&'^ — scrpix^^b^ 



and 



B' = ^ [e'^ + <&r + t^z^;. - - Vrbl^ [-apx'^bbrb^ + sdpiX^b%r + ^h6^] 

+ec0 [—(JpX^bbrbif, — 2(7pxx'bbrb^ — apx^b'brb^ — <Tpx^bb'j,b^ — apx^bbrb'^ 



+aiibl + 2aH&r&;] - 



sa'^x^b^br + 2s(j]ixx'b^br + 2saiiX^bb'br + saiiX^b^b'^. 



i-1 

c 



x'bbrb^bz + xb'brbijjbz + xbb'j.b^bz 
+xbbrb'^bz + xbbrb^b'^^ j [((Tq — txp) bb^_ + (Tpa;^6^ — sd'Ha;^^^^,/,] 
{—apx^bbrb 



- - 1 ) xbbrb^bz 



[{a'o-a'p)bbz 



bfj, + sa-ux^b^br + cr-ubl) 
+ (cto - ^p) b'bz + ((To - o-p) bb'z + o-pa;^&^ + 2CTpxa;'6^ + 'idpx^b'^b' ~ a^x'^b'^b^ - sa'^x'^ 
— 2s(7HXx'6^&0 — 2s(7iiX^bb'bcf, — saux^b^b'^]^ j [((tq — (Xp) bb^ + apx^fo^ — s(7Ha;^&^6</,] ^ 
Substituting Equation (|A5P into Equation 



which gives 



(|T7|) . we have 



- e^s] +eA - Ds 



l + Cs 
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TAYLOR SERIES EXPANSIONS FROM THE MIDPLANE 



As described in Section 12.1.31 the midplane is a singular point of the disk equations. The procedure adopted here 
for deahng with this issue is to calculate Taylor series expansions of the variables of integration from the midplane to 
a small distance As = 0.001 above the midplane before beginning numerical integration of the disk equations. The 
Taylor ser ies e xpansions for a general, non-ideal disk described using the conductivity tensor formalism are presented 
in Section fB.il The case of a disk in the ambipolar diffusion regime with constant neutral-field coupling parameter T 
is covered in Section [R2] 



General Case 

Differentiating Equation (|3T]) once with respect to s, we have x" — b'^x^ 

t" — h' 



3brX 



2^1 _ 



that vanish at the midplane gives 



b'j.x + 3bfx , dropping terms 



(Bl) 



Differentiating Equation (j48p once with respect to s and dropping terms that vanish at the midplane gives 

1" 



^2 - Vrfi - Vrfl) + 

+"^0,0 + {U,0 [K..0 + 1] - b'^,ojr,o) 



2 



(B2) 



-r.O ■ 



(^0 ^r.Q 



(B3) 



Differentiating Equation (j47p once with respect to s and dropping terms that vanish at the midplane gives 

'1 _ 5^ 

Differentiating the z component of Equation (15^1) gives e'^, q — e'^^Q+v'^f^b'^Q—Vrflb'^^Q = e^^nQ+e^^j^ q+w^ o^^r^o^^^o- 
Differentiating Equation (H^ yields e'^^j_ q = —b'^ QCcr.o — b'^ o^c0,o- Differentiating Equation (|55|) and Equation 



^P,oe'^z±ol /^o,a, respectively. These 



gives j;,o = (l/C-2)fe;^o and e^^n ,, = (j1,o + ^h.o 
relations can be used to solve for e^Q in Equation (jBSp . 

Calculating v'^q and v'J^ q is somewhat more complex. Differentiating Equations ((25)) and (P?]) twice with respect to 
s and dropping terms that vanish at the midplane gives 



+2 - Ox'^ + vx'^j^,^ - v-p'^j^fl + i/j^o - 2z/j;o^0 



(B4) 



and 



2 K,o 



«'-,0j ^0,0 = ^ ^ + 2W0,OW0,O - ^'./■.O {Vz,0 - Vrfi) 



Expressions for the unknown quantities j"o and j'^'q can be derived by differentiating Equations 
with respect to s and dropping terms that vanish at the midplane, which gives 



Jrfl 



and 



-f-O-p oGcnO + o'P,oecr±,o 



(B5) 

and ([55)1 twice 



(B6) 



We have 6q = b'^f^ + 6^0, since fe'^, q = 6"o = 0. Differentiating Equation P5)) twice gives the quantities e"^|| ^ and 
e"^ll Q as e"^|| = 26J, Qe[,^|| (, and e"^|| ^ = 26'^ Qe[,^||_o. The quantities e"^j^ q and e"^j^ q are obtained by differentiating 
Equation (|451) twice to obtain e"^j^ q = e"^ q ~ ^cril ^nd e"^^ q = e"^ q — e"^|| g. Finally, differentiating Equation ([39]) 
twice gives e"^ q = e"o + ^^^'^o ~ 2^2,0^0,0 and e"^ q = 2^2 0^0,0 — ^"0- Plugging these relations back into Equation (|B4I) 
and Equation (|B5[) gives two equations that can be solved for v'^q and v'^ q in terms of known quantities. 
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Ambipolar Diffusion Case with Constant T 
As in the general conductivity case, = 6J,q. 

Differentiating the relation B- {Vi — V) — once and evaluating at the midplane gives q = i 
b'ro {vir,o ~ ^r.o) + ^0 (''^i^.o ~ ^0,o)- Differentiating Equation (|63p twice and evaluating at the midplane yields 

•15^ 

Differentiating the z component of Equation (|8]) (cf. Equation (12) of L96) once and evaluating at the midplane 
gives 

/ fq™, _ 1 r oi \ 

(B9) 



')" - 1h' ?/ 

^ir,0 — ^"r.O^iz.O 



r,0 



(B8) 













[ 2 "-'H 







0<-l + ToK,o-<o])/^- 



Differentiating Equation (j64p three times and evaluating at the midplane yields 



" "<Pfi II t I Ul 1.1 \ 



(BIO) 



"r,0 



Finally, differentiating the r and </> components of Equation ([5]) (cf. Equations (13) and (14) of L96) twice and 
evaluating at the midplane gives 



II I O -v^ II I -»^// / // " \ 



and 



These equations can be solved for w"o and v'^ q in terms of known quantities. 



2 



(Bll) 



(B12) 



18 



REFERENCES 



Andrews, S. M., & Williams, J. P. 2007, ApJ, 659, 705 

Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 

Bally, J., Reipurth, B., & Davis, C. J. 2007, in Protostars & 

Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: 

Univ. Arizona Press), 215 
Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 (BP82) 
Brauer, F., DuUemond, C. P., & Henning, Th. 2008, A&A 480, 

859 

Campbell, C. G. 2003, MNRAS, 345, 123 
Campbell, C. G. 2005, MNRAS, 361, 396 
Campbell, C. G. 2009, MNRAS. 392, 271 
Cao, X., & Spruit, H. C. 2002, A&A, 385, 289 
Casse, F., & Ferreira, J. 2000, A&A, 361, 1178 
Casse, F., & Keppens, R. 2002, ApJ, 581, 988 
Casse, F., & Keppens, R. 2004, ApJ, 601, 90 
Ciolek, G. E., & Mouschovias, T. Ch. 1993, ApJ, 418, 774 
Draine, B. T., Roberge, W. G., & Dalgarno, A. 1983, ApJ, 264, 
485 

Ferreira, J. 1997, A&A, 319, 340 (F97) 

Ferreira, J., & Casse, F. 2004, ApJ, 601, L139 

Guan, X., & Gammie, C. F. 2009, ApJ, 697, 1901 

Hartigan, P., Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736 

Hayashi, C. 1981, Prog Theor Phys Supp, 70, 35 

Heyvaerts, J., & Norman. C. 1989, ApJ, 347, 1055 

Keller, Ch., & Gail, H.-P. 2004, A&A. 415. 1177 

Konigl, A. 1989, ApJ, 342, 208 

Konigl, A. 2004, ApJ, 617, 1267 (K04) 

Konigl, A., & Salmeron, R. 2011, in Physical Processes in 
Circumstellar Disks around Young Stars, ed. P. J. V. Garcia 
(Chicago: Univ. Chicago Press), in press (arXiv: 1004: 1875) 



Konigl, A., Salmeron, R., & Wardle, M. 2010, MNRAS, 401, 479 

Lesur, G., & Longaretti, P.-Y. 2009, A&A, 504, 309 

Levy, E. H., & Sonett, C. P. 1978, in Protostars & Planets, ed. T. 

Gehrels (Tucson: Univ. Arizona Press), 516 
Li, Z.-Y. 1995, ApJ, 444, 848 
Li, Z.-Y. 1996, ApJ, 465, 855 (L96) 

Lovelace, R. V. E., Rothstein, D. M., & Bisnovatyi-Kogan, G. S. 

2009, ApJ, 701, 885 
Lubow, S. H., Papaloizou, J. C. B., & Pringle J. E. 1994, 

MNRAS, 267, 235 
Lubow, S. H., Papaloizou, J. C. B., & Pringle J. E. 1994, 

MNRAS, 268, 1010 
Nishi, R., Nakano, T., & Umebayashi, T. 1991, ApJ, 368, 181 
Ogilvie, G. L & Livio, M. 2001, ApJ, 553, 158 (OLOl) 
Guyed, R., & Pudritz, R. E. 1997, ApJ, 482, 712 
Rothstein, D. M., & Lovelace, R. V. E. 2008, ApJ, 677, 1221 
Salmeron, R., Konigl, A., & Wardle, M. 2007, MNRAS, 375, 177 

(SKW07) 

Salmeron, R., Konigl, A., & Wardle, M. 2011, MNRAS, in press 

(arXiv:1006.0299) (SKWll) 
Salmeron, R., & Wardle, M. 2005, MNRAS, 361, 45 
Shu, F. H., Shang, H., Gounelle, M., Glassgold, A. E., & Lee, T. 

2001, ApJ, 548, 1029 
Vlahakis, N., Tsinganos, K., Sauty, C, & Trussoni, E. 2000, 

MNRAS, 318, 417 
Wardle, M. 2007, Ap&SS, 311, 35 

Wardle, M., & Konigl, A. 1993, ApJ, 410, 218 (WK93) 
Zanni, C, Ferrari, A., Rosner, R., Bodo, G., & Massaglia, S. 
2007, A&A, 469, 811 



